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SECTION  A.  INTRODUCTION 


This  report  is  concerned  with  the  mathematical  simulation  of 
the  motion  of  an  occupant  of  a  motor  vehicle  which  has  been  involved 
in  an  accident.  There  are  three  phases  to  this  simulation  project. 

The  first  phase  is  quite  similar  to  the  Cornell  simulation  project, 
differing  only  in  the  methods  of  obtaining  certain  quantities  and  in 
the  number  of  allowed  body  contacts  with  the  vehicle  interior.  An 
Highway  Safety  Research  Institute  computer  program  is  currently  being 
used  to  generate  results  which  are  being  compared  with  results  from  the 
Cornell  program. 

The  second  phase  is  the  main  concern  of  this  report.  This  new 
model  includes  all  changes  and  improvements  of  the  earlier  model  which 
can  be  made  without  increasing  the  number  of  degrees  of  freedom  of  the 
system.  The  vehicle  is  equipped  with  a  restraint  system  as  shown  in 
Fig.  1.  The  motion  considered  occurs  only  in  the  forward  vertical  plane 
(i.e.  in  two  dimensions)  as  a  result  of  either  forward  acceleration  (say 
from  a  rear-end  collision)  or  deceleration  from  a  front-end  collision 
or  sudden  stop.  A  rear-seat  passenger  and  injury  potential  are  considered 
in  this  model.  It  has  not  yet  been  programmed  for  the  digital  computer 
but  algorithms  for  doing  so  are  laid  out  in  this  report. 

The  third  phase  of  the  simulation  project  proceeds  in  two  directions. 
The  first  is  the  development  of  simpler  models  for  the  study  of  specific 
aspects  of  the  collision  problem.  The  second  proceeds  in  the  opposite 
direction  with  the  objective  of  further  sophisticating  the  model  by  the 
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Figure  1 

THE  SYSTEM  -  FRONT  SEAT 
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addition  of  body  elasticity,  three-dimensional  motion,  asymetries,  and 
relative  vehicle  interior  motions. 

The  body  of  the  person  is  approximated  by  an  eight  segment 
articulated  system  of  links.  The  eight  angles  of  these  segment  center 
lines,  which  are  measured  counter-clockwise  from  the  true  horizontal, 
and  the  two  inertial  Cartesian  coordinates  of  the  hip  joint  comprise  the 
ten  degrees  of  freedom  of  the  body  in  the  simulations  of  phases  one 
and  two.  Phase  three  will  have  a  larger  number  of  degrees  of  freedom  due 
to  the  explicit  simulation  of  such  body  deformabilities  as  spine-limb 
length,  shoulder  joint  motion,  and  chest- face-abdomen  deformations.  The 
phase  one  program  and  Cornell  ignore  these  deformabilities  while  the 
phase  two  simulation  approximates  them  by  the  introduction  of  injury 
criteria  based  on  forces.  The  chest-face-abdomen  deformabilities  are 
further  approximated  by  the  use  of  force-deflection  loading  curves  which 
modify  their  contact  or  belt  forces.  See  Section  C.4.e  and  C.4.f  for  details. 
These  injury  criteria  can  also  be  used  as  end-of-program  tests,  with  a 
program  option  to  specify  which  of  these  criteria  would  be  so  used.  These 
will  be  enumerated  below.  In  phases  two  and  three,  each  of  the  seven  joints 
of  the  body  is  considered  to  possess  torsional  elasticity  proportional  to 
the  relative  joint  angle,  Coulomb  friction,  and  2-sided  angular  restraints 
approximating  actual  limitations  to  joint  motion.  Cornell  allows  elasticity 
for  only  the  two  spinal  joints,  with  two-sided  linear  restraints  for  the 
neck  and  spinal  joints  and  one-sided  linear  restraints  for  the  elbow 
and  knee.  The  phase  one  program  adds  a  two-sided  linear  restraint 


for  the  shoulder  and  a  one-sided  one  for  the  hip.  The  phase  two 
simulation  allows  the  form  of  the  restraints  to  be  linear,  cubic,  or 
arbitrary.  Again  an  option  in  the  computer  program  will  specify  which 
of  these  would  be  used.  See  Section  C.A.c  for  details. 

The  vehicle  of  phase  one  and  the  Cornell  vehicle  possess  one  degree 
of  freedom.  This  is  straight-line  horizontal  motion.  Vertical  motion 
and  pitch  are  included  in  the  models  of  phases  two  and  three.  In  the 
former  simulation  an  option  exists  which  specifies  whether  this  horizontal 
acceleration  has  an  arbitrarily  specified  time  history  or  whether  a  cart- 
stopping-mechanism  is  used  as  would  be  the  case  with  a  sled  test.  The 
phase  two  simulation  has  no  such  option  and  all  three  accelerations  must 
be  known  time-functions.  However,  it  is  possible  to  include  tables  for 
the  horizontal,  vertical,  and  angular  positions  and/or  rates  as  input  data 
and  thus  eliminate  some  or  all  of  the  vehicle  integrations.  The  phase 
three  simulation  will  have  the  option  of  using  explicit  differential 
equations  for  the  three  vehicle  variables,  a  three-dimensional  cart-stopping- 
mechanism,  or  the  above  tabular  data. 

The  vehicle  interior  is  described  by  taking  several  vertical  sections 
through  the  vehicle,  through  the  center  of  the  occupant,  through  his  legs, 
and  through  his  arms.  These  sections  are  then  superimposed  to  form  one 
composite  interior  layout,  as  shown  in  figure  2.  All  three  phases  have 
an  option  to  specify  the  position  of  the  person  in  the  vehicle,  whether 
driver,  front-seat  passenger  or  rear-seat  passenger.  Cornell  specifies 
only  driver  or  front-seat  passenger  by  the  values  of  certain  contact 
surface  parameters.  The  Cornell  vehicle  interior  contains  only  seven 
contact  surfaces,  all  of  which  are  flat.  The  phase  one  program  has  twelve 
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flat  surfaces ,  while  phases  two  and  three  have  seventeen,  some  of  which 
may  be  curved,  Cornell  considers  only  seven  combinations  of  contact 
between  particular  body  segments  and  contact  surfaces.  The  number  of 
possible  contacts  in  the  Highway  Safety  Research  Institute  simulations 
varies  with  the  position  of  the  occupant  in  the  car.  In  the  phase  one 
program  the  driver  has  17  possible  contacts,  the  front-seat  passenger  14, 
and  the  rear-seat  passenger  13.  In  phases  two  and  three  these  numbers 
are  20,  18  and  14  respectively. 

The  Cornell  contact  surfaces  are  assumed  to  consist  of  a  single 
material.  This  is  also  the  case  in  the  phase  one  program.  In  phases 
two  and  three  some  of  the  contact  surfaces  may  be  padded,  as  indicated 
In  Table  10.  Cornell  and  the  phase  one  program  have  only  rigid  (non¬ 
movable)  contact  surfaces,  while  in  phase  two  the  motion  of  the  non-rigid 
ones  (for  example,  the  sunvisor,  steering  wheel  rim  and  column,  and  head 
rest)  is  approximated  by  force-deflection  loading  curves  and  the  use  of 
the  shared  deflection  subroutine  (Fig.  23).  In  phase  three  the  non-rigid 
motion  will  be  expressed  in  terms  of  explicit  differential  equations.  All 
these  force-deflection  loading  curves  are  assumed  to  be  expressed  by 
fifth-order  polynomials  in  both  the  deflection  and  the  rate  of  deflection. 

Cornell  and  the  phase  one  program  represent  the  seat  by  two  systems, 
one  a  linear  spring  whose  force  acts  at  the  front  edge  of  the  seat  and  the 
other  a  linear  damper  -  third  order  non-linear  spring  combination  whose 
force  acts  at  the  hip  joint.  Thus  this  latter  force  may  move  along  the 
seat  following  body  motion  and  might  overlap  the  front  edge  force.  This 
overlap  Is  prevented  in  phases  two  and  three  either  by  logic  or  by  using 
an  entirely  different  seat  representation,  two  possibilities  for  which  are  a 
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rigid  beam  supported  by  the  present  spring-damper  systems  and  a  pair 
of  non-rigid  beams.  In  all  cases,  logic  allows  the  person  to  fall 
completely  off  the  seat.  Logic  in  the  appropriate  sections  of  the 
simulations  of  all  three  phases  automatically  differentiates  between 
a  person  whose  feet  touch  the  floor  and  a  child  whose  legs  stick  straight 
out  in  front  of  him. 

In  all  cases,  the  restraint  system  is  composed  of  a  lap  belt  whose 
force  acts  at  a  constant  radial  distance  from  the  hip  joint  and  a  shoulder 
restraint  whose  two  forces  act  at  specified  points  on  the  body  torso 
center  lines.  Cornell  has  an  option  to  specify  whether  the  3  belt  sections 
and  thus  the  three  belt  forces  are  independent,  whether  the  lap  belt  is 
independent  of  the  shoulder  restraint  while  the  tension  forces  in  the 
two  sections  of  the  shoulder  restraint  are  equal,  or  whether  the  forces 
in  all  three  sections  are  the  same.  The  phase  one  program  and  the  phase 
two  simulation  assume  that  the  first  of  these  alternatives  is  true.  Phase 
three  will  have  this  option.  All  three  phases  have  another  option  to 
specify  how  many  of  the  belts  are  in  use,  whereas  Cornell  must  use  the 
values  of  the  various  belt  parameters  to  do  this.  The  phase  one  program 
assumes  that  the  chest  wall  and  abdomen  do  not  deform,  while  phase  two 
approximates  their  deformation  by  load-deflection  curves  and  a  shared- 
deflection  subroutine  to  generate  the  belt  forces.  Phase  three  will  have 
explicit  differential  equations  for  their  motion. 

Cornell  and  the  phase  one  program  do  not  consider  injuries  at  all. 

In  phase  two,  the  injuries  which  are  considered  in  the  first  approximation 
for  non-rigid  body  segments  are  joint  dislocations,  bone  fractures,  and 
cuts.  Joint  dislocations  can  be  of  two  types;  one  caused  by  excessive 
relative  angle  and  the  other  by  excessive  tensile  forces  across  the  joint. 
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Fractures  are  also  of  two  types;  one  caused  by  excessive  compressive 
force  along  the  center  line  of  the  bone,  the  other  by  excessive  shear 
forces.  This  last  type  affects  not  only  the  spine  and  limbs,  but  also 
the  skull,  face,  and  ribs.  Cuts  are  considered  only  for  head-windshield 
contact  when  the  glass  is  actually  shattered.  These  injury  criteria 
will  also  appear  in  phase  three,  though  perhaps  in  different  forms. 

For  all  three  phases  the  equations  of  motion  for  the  system  are 
obtained  from  the  system  Lagranglan,  which  includes  the  system  kinetic 
energy,  potential  energy,  dissipated  energy  rate,  and  the  classical 
generalized  non-conservative  forces. 


SECTION  B. 


Derivation  of  Matrix  Equation 

00 

The  body  matrix  equation,  Az  =  b,  where  z  is  the  generalized  co¬ 
ordinate  vector  with  the  ten  components  through  ag,  x^  and  y^;  A  is 
the  body  matrix  as  shown  in  section  C.5;  and  b  is  the  generalized  force 
vector  whose  components  are  calculated  in  section  C.4,  is  obtained  from 
the  system  Lagrangian: 


<L  ,  3(KE)  _  KKEl  +  iM  +  1M  =  F  ,  for  i  =  1-10, 

dt  L  J  zi 

where  KE  is  the  system  kinetic  energy,  PE  is  its  potential  energy,  DE 
is  its  dissipated  energy  rate,  and  the  are  the  classical  generalized 
forces,  developed  in  detail  in  sections  C.4.e  and  C.4.f. 

The  total  body  kinetic  energy  is 


where  m  and  I.  are  the  mass  and  moment  of  inertia  of  the  i  th  body  se- 
i  1 

gment ,  respectively,  cl  is  the  angle  of  the  i-th  body  segment  center 

line  to  the  true  or  inertial  horizontal  measured  counter-clockwise,  and 

X  Y  are  the  inertial  coordinates  of  the  center  of  gravity  of  the  i-th 
i’  i 

body  segment  as  follows: 


1  "  *h  +  Pl 


+  p.,  cosa, 


Y1  =  yh  +  pi  sinax 


xi  =  *h  “  piai  sin“i 

^i  =  yh  +  cosa1 


9 


X2  «  xh  +  cosc^  +  P2  cosa2  ^2  =  *h  ~  Ll“l  sinal  "  p2a2  8ina2 

Y2  =  yh  +  \  sinc^  +  P2  slnct2  *2  “  K  +  Ll“l  C°Sal  +  P2°2  COSa2 

X3  «  ^  cosa^  +  L2  cosa2  +  P3  cosa3 

Y3  =  yh  +  L.^  sina^  +  L2  sina2  +  P3  sina3 

X3  -  sina^  -  L2a2  sina2  -  P3a3  slna3 

Y3  =  yh  +  cosc^  +  L2a2  cosa2  +  p^  cosa3 

X4  -  xh  +  Lx  cosc^  +  L2  cosa2  +  L3  cosa3  +  P4  cosa4 
Y4  =  yh  +  L3  sina1  +  L2  slna2  +  L3  sina3  +  P4  sina4 

X4  =  x^  -L^a^  sina^  _L20t2  sin0t2  -L3a3  sina3  _p4a4  sina4 
Y4  =  yh  +  cosa^  +  L2<x2  cosa2  +  L3a3  cosa3  +  P4a4  cosa4 

X^  =  x^  +  cosa^  +  L2  coso2  +  L4  cosa3  +  p^  cosa^ 

Y5  =  yh  +  Lx  sina1  +  L2  sina2  +  sina3  +  P5  sina5 

X,.  =  x^  -L^a^  sinai  -L2a2  sina2  ~L4a3  sina3  -p5a5  sina5 
Y5  =  y^  +  cosa1  +  L2a2  cosa2  +  I^a-j  cosa3  +  P5a5  cosa5 

X6  =  xh  +  \  cosa1  +  L2  cosa2  +  cosa3  +  L5  cosa5  +  Pfe  cosa6 

Y6  =  yh  +  Lx  sina1  +  L2  sina2  +  L4  sina3  +  sina5  +  p&  sin«6 

X5  =  Xn  -L.^  sina1  -L^  sina2  -L^  sina3  +  L5a5  sina5  -P6a6  slnct6 
Y6  =  yh  +  L.^  cosa1  +  cosot2  +  cosa3  +  L5a5  cosa5  +P6“6  cosa6 
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+  Py  cosa^  X 

yh  +  P7  sina7  Y 

x,  +  L-,  cosa-  +  p0  cosa0 

a  /  /  o  o 

x,  +  L  sina-  +  pQ  sinaQ 
n  /  /  o  o 


xh 

-p^a^  sina^ 

• 

yh 

• 

+p^a7  cosa^ 

• 

•  • 

xn 

=  x,  -  L-,a., 

sina_, 

8 

h  7  7 

7 

• 

•  • 

00 

“  yh  +  L?a? 

cosa^ 

-p8a8  sinag 

+  Pg^g  cosag 


Figures  3,4  and  5  show  the  body  dimensional  parameters  and  angles;  is 
the  link  length,  p^  the  distance  from  the  previous  joint  to  the  center  of 
gravity,  and  as  above.  In  figure  3  the  axes  shown  with  broken  lines 
are  the  true,  outside,  or  inertial  horizontal  and  vertical  references, 
while  the  solid  line  axes  are  relative  to  the  interior  of  the  vehicle  and 
are  thus  inclined  at  an  angle  $  to  the  inertial  axes.  This  angle  $  is 
the  vehicle  pitch  angle  and  is  positive  for  a  vehicle  having  its  front 
end  down.  The  angles  of  the  body  segment  center  lines  with  the  vehicle 
internal  horizontal  ref erence(i.e.  the  floor)  are  equal  to  the  inertial 
angles  plus  the  vehicle  pitch  angle;  i.e.  0^  -  a^  +  $ .  If  the  vehicle 
does  not  pitch,  then  the  reference  frames  will  coincide.  In  all  cases, 
x^  and  y^  are  the  inertial  coordinates  of  the  hip  joint. 

The  potential  energy  due  to  joint  elasticity,  the  seat  springs  (both 
linear  and  third-order  non-linear)  and  gravity,  respectively,  is: 


PE 


i=l 


4k.  (a. -a  +a  -a.  ) ^  +  4s(y1+y  -y  )^ 
2  i  l  m  mo  io'  2  J1  Jz  J  zo' 
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Figure  3 

BODY  ANGLES  AND  PITCH 
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Figure  5 

BODY  CONTACT  DIMENSIONS 


where  y^  and  y^  are  auxilliary  variables  to  be  defined  in  section  C.4.d. 
In  general,  potential  energy  is  the  integral  of  the  force  or  torque 
over  the  distance  (linear  or  angular)  involved.  Since  these  are  all 
well  defined  functions  in  the  analysis,  they  can  all  be  directly  in¬ 
tegrated;  even  the  seat  non-linear  spring,  which  gives 


3 

I 


e.y 


j+i 

l 


J+i 


However,  if  either  of  the  other  two  seat  representations  are  used, 
appropriate  terms  must  be  substituted. 

The  dissipated  energy  rate  due  to  seat  friction,  the  seat  damper, 
joint  friction,  and  the  joint  restraints,  respectively,  is: 


DE  -  /  f 
0 


7  <t+a 


T  .da 
i  m 


l  T^i] 


In  each  case  here,  the  rate  of  energy  dissipation  is  equal  to  the  integral 
of  the  force  or  torque  over  the  velocity  (linear  or  angular) .  In 
particular  if  the  joint  restraints  (section  C.4.c.  and  Fig.  9)  do  not 
depend  on  the  angular  rates,  they  could  be  integrated  yielding: 


ZVvV- 
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The  F  .  are  obtained  by  finding  the  external  forces  and  moments 
zi 

consisting  of  the  three  belt  forces,  the  normal  and  tangential  contact 
forces,  and  the  pitch.  (Sections  C.4.e.  and  C.4.f.). 

After  differentiation  and  proper  algebraic  manipulation,  there 
results  a  set  of  ten  differential  equations  in  the  ten  generalized 
coordinates,  which  can  be  recognized  as  A  z  =  b,  where  A  is  a  function  of 
z  and  b  of  z  and  z.  For  ease  of  calculation  b  is  repartitioned  into 
those  items  calculated  in  Section  C.4. 
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SECTION  C. 


SIMULATION 


The  detailed  simulation  presented  in  this  section  is  that  of 
phase  two.  It  incorporates  all  the  additions  and  changes  that  can 
be  made  to  the  phase  one  program  without  raising  the  number  of  degrees 
of  freedom  of  the  system.  The  overall  flow  chart  for  the  phase  two 
simulation  program  is  shown  in  figure  6. 


SECTION  C.l 
INPUT  INFORMATION 

The  general  format  for  the  input  information  is  that  the  various 

parameters  should  be  grouped  according  to  type.  This  arrangement  allow 

any  one  section  to  be  changed  by  simply  substituting  new  data  cards. 

A  complete  list  of  parameters  and  the  basic  variables  along  with  their 

units  and  definitions  will  be  found  just  before  the  introduction,  section  A. 

Arbitrarily  specifiable  initial  conditons  are:  body  position 

angle  (0,  through  0  ),  belt  length  (i  )  and  angle  (6)  for  j'«  1 

lo  8o 

through  3,  vehicle  velocity  (x  ). 

co 

Vehicle  pitch  angle,  plus  its  rate  and  acceleration,  and  deceleration 
time-histories  are  entered  in  tabular  form,  as  well  as  vehicle  velocity 
and  position,  if  available.  [xc(0,  yc(0,  with  xc(tm)  *  yc(tm^> 

$(tm),  xc(tm),  yc(tm),  i(tm)  if  available,  for  all  tm]  The  only  other 
vehicle  parameters  needed  at  present  are  its  mass  (mc)  and  moment  of 
inertia  (Ic). 
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Figure  6 

OVERALL  FLOW  DIAGRAM 
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subscript 

i 

n 

2 

3 

4 

5 

6 

7 

8 

9  10 

other 

contact  radius 

r 

X 

X 

X 

X 

X 

X 

X 

X 

X  X 

rA 

shear  fracture 
force 

FSF 

X 

X 

X 

X 

X 

X 

X 

X 

X 

mass 

m 

X 

X 

X 

X 

X 

X 

inertia 

I 

X 

X 

X 

X 

X 

X 

length 

L 

X 

X 

X 

X 

X 

X 

X 

X 

center  of  gravity 
distance 

P 

X 

X 

X 

X 

X 

X 

X 

X 

tension  dislocation 
force 

FTD 

X 

X 

X 

X 

X 

X 

X 

compression  fracture 
force 

FCF 

X 

X 

X 

X 

X 

X 

X 

coulomb  friction 
coefficient 

C* 

X 

X 

X 

X 

X 

X 

X 

elasticity 

K 

X 

X 

X 

X 

X 

X 

X 

friction  velocity 
limit 

e 

X 

X 

X 

X 

X 

X 

X 

minimum  joint  angle 

U) 

X 

X 

X 

X 

X 

X 

maximum  joint  angle 

ft 

X 

X 

X 

X 

X 

X 

linear  joint  stop>  f 

T* 

X 

X 

X 

X 

X 

X 

coefficients  j 

Tf 

X 

X 

X 

X 

X 

X 

center  of  curvature 
distance 

p1 

X 

X 

X 

center  of  curvature 
angle  from  center- 
line 

Y 

1 

X 

TABLE  1. 
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Parameters  for  the  body  (subscript  i),  belts  (subscript  j), 

contact  surface  positions  (subscript  k)  and  materials  (subscript  i) 

are  shown  in  tables  1  through  4,  respectively.  Table  5  shows  the 

various  meanings  attached  to  the  body  subscript  i  for  the  body 

parameters.  The  present  seat  parameters  are:  cg>  s,  $3*  zo’ 

z,  y  ,  W  .  and  F1  . 

*  fO  o’  so 


subscript  j 


parameter 

1.  lap  belt 

2.  lower  torso 

3.  upper  torso 

h 

X 

X 

X 

l 

X 

X 

X 

0 

4> 

X 

X 

X 

0 

V 

X 

0 

X 

0 

A 

X 

X 

X 

N 

X 

X 

X 

G 

X 

X 

X 

R 

X 

X 

X 

0 

X 

X 

X 

n 

°n+N 

? 

? 

? 

TABLE  2. 


subscript 

k 

x" 

yM 

parameters 

V  D  D"  other 

1 

head  rest 

X 

X 

X 

X 

X 

2 

windshield 

X 

X 

X 

X 

X 

3 

sun  visor 

X 

X 

X 

X 

X 

4 

upper  panel 

X 

X 

X 

5 

lower  panel 

X 

X 

X 

X 

X 

6 

toeboard 

X 

X 

X 

X 

X 

7 

back  of  seat 

X 

X 

X 

X 

X 

8 

column 

X 

X 

X 

X 

X 

9 

seat  back 

X 

X 

X 

X 

*  Dfr 

10 

rim 

X 

X 

X 

x  Dt,Ds 

11 

roof 

X 

12 

floor 

X 

TABLE  3. 
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subscript 


parameters 


i 

M 

e 

G 

R 

N 

0 

n 

°n+N 

other 

1 

•eat 

X 

X 

X 

X 

X 

X 

? 

2 

glass 

X 

X 

X 

X 

X 

X 

? 

FSW 

3 

padding 

X 

X 

X 

X 

X 

X 

? 

n 

4 

floormat 

X 

X 

X 

X 

X 

X 

? 

n 

5 

floor 

X 

X 

X 

X 

X 

X 

? 

6 

metal 

X 

X 

X 

X 

? 

7 

rim 

X 

X 

X 

X 

X 

X 

? 

8 

headrest 

X 

X 

X 

X 

X 

X 

? 

9 

visor 

X 

X 

X 

X 

? 

10 

column 

X 

X 

X 

X 

? 

11 

abdomen 

X 

X 

X 

X 

? 

12 

chest 

X 

X 

X 

X 

? 

13 

face 

X 

X 

X 

X 

? 

14 

shoulder 

X 

X 

X 

X 

7 

TABLE  4. 
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subscript 

bone  injury 

criteria 

i 

body 

segment 

loint 

contact 

arc 

shear 

fracture 

compression 

fracture 

parameters 

?\K,T\ 

t  , 

E,  FTD 

*.p\y 

FSF 

FCF 

1 

lower 

torso 

hip 

lower 

back 

lower 

spine 

lower 

spine 

2 

center 

torso 

lower 

spine 

upper 

back 

center 

spine 

center 

spine 

3 

upper 

torso 

upper 

spine 

chest 

upper 

spine 

upper 

spine 

4 

head 

neck 

head 

skull 

5 

upper 

arm 

shoulder 

shoulder 

back 

upper 

arm 

6 

lower 

arm 

elbow 

elbow 

rib 

lower 

arm 

7 

upper 

leg 

knee 

knee 

upper 

leg 

upper 

leg 

8 

lower 

leg 

foot 

lower 

leg 

lower 

leg 

9 

shin 

face 

10 

i 

hand 

TABLE  5. 
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SECTION  C.2 


PRELIMINARY  CALCULATIONS 


In  the  program  for  the  simulation,  certain  parameters  are  computed 
here  and  stored  for  later  use,  so  as  to  keep  the  amount  of  calculation 
to  a  minimum.  One  set  comprises  the  parameters  used  in  calculating 
the  elements  of  the  body  matrix  (a^  fcr  i  =  1  to  17);  another  the 
initial  belt  parameters  (X^  for  i  ■  1  to  6) .  Their  defining  relations 


are  as  follows: 


i  =  4,6,8 
i  =  5,7 

i  =  1,2 


i  =  10,11,12 


i  =  13,15,17 
i  =  14,16 


a.  =  m.p. 
i  i  i 


a.  =  Vi  +  mi+iLi 


a  =  m  p  +  L  l  m 
1  11  j=i+l  J 

a3  =  m3P3  +  mAL3  +  (m5+m6)LA 


8 

aQ  =  l  m. 

9  J-l  3 

2  .2  § 

ai  =  1i-9  +  mi-9  p i-9  +  Li-9j£._8  . 

2  2  2 
a12  =  X3  +  3  +  m4L3  +  (m5+m6)L4 


m 


ai  =  Ii-*9  +  mi-9pi-9 


a-  =  1-*  9  +  mi  9^i  9  ^  Li-9mi-8 


i  l- 
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X1  ■  lo  c08<t’o  *  A2  “  lo  ,in*o 


x3  “  L18in®lo  +  h2sln02o  "  Sosin<fl2o 


LlCos0lo  +  .h2cos02o  -  ^2ocos^o 


Xc  *  L^sinQ-  +  Losin0„  +  hosin0o  -  sin<f)0 
5  1  lo  2  2o  3  3o  3o  T3o 


L-cos0-  +  Locos0«  +  tucos©,, 

1  lo  2  2o  3  3o 


So  COSSo 


Single  parameters  are  the  initial  lap  belt  angle  and  length,  and  the 
angle  between  the  shin  surface  and  the  lower  leg  center  line.  These 
are  defined,  respectively,  as: 

ho  =  *0  +  sin_1[hi/Ao]' 

So  -  S2-hl  +  hl[*o  -  0lo  +  2  n].’ 

%  =  sin_1[(r7-rA)/L8]. 

Other  calculations  performed  here  are  of  the  initial  values  of 
the  output  data  variables.  Most  of  these  will  be  zero;  eg  the  head 
and  chest  accelerations  at  their  centers  of  gravity,  all  energies 
except  kinetic,  and  all  contact  forces  except  those  on  on  the  seat 
back,  head  rest,  floor,  and  toeboard.  Initial  kinetic  energies  are, 
respectively  for  the  vehicle  and  body  segments: 


KE 


rm  x  , 
2  c  co’ 


KE.  =  Thn.  x^  ,  for  i 
i  2  i  co  * 


1-8. 


Initial  contact  forces  are  dependent  on  initial  body  position  (or 
perhaps  vice  versa, using  the  appropriate  loading  relation  for  force 
in  terms  of  deflection  with  rate  of  deflection  equal  to  zero). 
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Lastly  the  initial  accelerations  on  the  body  are  calculated 
considering  only  the  initial  forces  due  to  gravity,  the  belts,  and 
the  initial  contacts.  The  generalized  centrifugal  forces,  joint 
torques,  and  net  seat  cushion  forces  are  all  zero. 

The  vehicle  initial  accelerations  are  read  from  the  tables. 
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SECTION  C.3 


INTEGRATION 


At  present,  a  four  point  Adams -Moult on  method  started  with  a 
modified  Euler  method  is  used.  The  Euler  method  is  used  to  obtain 
the  first  three  points  after  the  given  initial  point.  Thereafter 
four  points  are  used  in  the  Adams  third  difference  prediction  and 
correction  equations.  At  each  computation  of  a  new  point,  the 
iteration  continues  until  successive  approximations  are  within  a 
given  tolerance  or  the  iteration  has  been  tried  a  prescribed  number 
of  times. 

With  an  integration  option,  different  schemes  could  be  chosen 
according  to  the  option.  Otherwise  the  integration  would  be  (and  is) 
carried  out  by  a  subroutine. 

If,  in  the  input  information,  tabular  data  is  available  only 

for  Xcftm)s  yc(tm>>  and  then  these  must  be  inteErated  for  the 

vehicle  position  and  velocity  and  for  its  pitch  angle  and  rate.  If 

tabular  data  is  also  available  for  xc(tm),  yc^m^’  ^fcm^  and/°r  for 

x  (t  ),  v  (t  ),  *(t  ),  some  or  all  of  the  integrations  for  the 
c  m  *  J  c  m 9  nr 

vehicle  may  be  omitted  and  the  proper  values  read  in  from  the  table. 

The  body  generalized  coordinates  (x^,  yh>  -  ag)  are  now  in  the 
inertial  reference  frame  and  can  be  directly  integrated  twice.  The 
body  angles  and  rates  relative  to  the  vehicle  interior  [e^t+At), 


8^(t+At)]  are  obtained  by  adding  *(tnrfl)  t0  «(t+At)  and 

a(t+At)  respectively. 

Figure  7  shows  the  various  reference  frames,  along  with  the 
vehicle  center  of  gravity  and  body  positions  at  some  later  time  t»tm* 


28 


REFERENCE  FRAMES 


SECTION  C.4 


GENERALIZED  FORCES 

Once  the  set  of  ten  equations  in  the  ten  generalized  coordinates 

is  obtained  for  the  body  motion  (see  Section  B  for  details),  it  can 

be  separated  so  that  all  the  acceleration  terms  are  on  one  side;  i.e. 

**  • 

Az  *  b,  where  A  is  a  function  of  z  and  b  of  z  and  z.  In  order  to 
facilitate  the  calculation  of  the  components  of  the  vector  b,  which  is 
now  called  the  generalized  force  vector,  it  is  partitioned  into  sections 
for  the  forces  and  moments  due  to:  1)  centrifugal  force ,  2)  gravity, 

3)  joint  torques,  4)  the  seat  cushion,  5)  the  belts,  and  6)  contacts 
with  the  interior  surfaces  of  the  vehicle,  as  shown  in  Fig.  8.  In 
order  then, 

b^  *  Cj-  +  Q^,  for  i  *  1  to  10. 

The  vector  s  has  as  components  through  a^,  x^,  and  y^. 

In  the  present  case,  when  the  magnitude  of  the  vehicle  deceleration 
is  known  for  each  time  point,  the  magnitude  of  the  resulting  force  on 
the  vehicle  is  simply: 

*c(tm)  -  ♦  yin*)' 

The  instantaneous  positions  and  velocities  of  the  seven  joints, 
relative  to  the  inside  of  the  vehicle,  are  needed  in  several  of  the 
subsequent  sections.  They  are  calculated  from  the  following  relations. 
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31 


hip 


i=l  X-^t)  =  [xh(t)-xc(t) ]cos$(tm)-[yh(t)-yc(t) ]sin$(tm) 

yi<t)  =  [xh(t)-xc(t)]sin$(t]n)  +  [yh(t)-yc(t)]cos4>(tm) 

and 

x  (t)  =  [x,  (t)-x  (t)]cos$(t  )-[y  (t)-y  (t)]sin$(t  )-y  (t)i(t  ) 
1  nc  mnc  ml  m 

y.d:)  =  [x  (t)-x  (t)]sin$(t  )+[y  (t)-y  (t)]cos$(t  )+x  (t)$(t  ) 
1  nc  m  n  c  ml  m 

knee 

i=7  x^(t)  =  x^(t)+L^cos0^ (t) 

y7(t)  =  y1(t)+L7sin07(t) 

x? d)  =  x1(t)-L707(t)sin07(t) 
y7(t)  =  y1(t)+L707(t)cos07(t) 
torso,  neck,  elbow 

i=2,3  x1(t)  =  x±_1(t)  +  Li_1cos0i_1(t) 

^’6  y±d)  =  y^d)  +  Li_1sin0il(t) 

Id)  =  i^d)  -  ^(tisinO^jCt) 

yAd)  =  y1_1(t)  +  L1_10i_1(t)cos01_1(t) 

shoulder 

i=5  x,.(t)  “  +  L^cosG^(t) 

y5(t)  =  y3(t)  +  L^sinO^Ct) 

x5(t)  =  x3(t)  -  L^03(t)sin03(t) 
y5(t)  =  y3(t)  +  L^03(t)cos03(t) 

where  0. (t)  =  a.(t)  +  $(t  ). 

11  m 
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SECTION  C.4.a 
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SECTION  C.4.b 


GRAVITY 

Torque  due  to  gravity  in  the  direction  is: 

G±(t)  =  gaicosai 

for  i  =  1  to  8,  where  the  a^  were  defined  and  calculated  in  section 
C.2.  In  the  true  horizontal  the  force  is 
Gg(t)  =  0; 

in  the  true  vertical  direction,  it  is 

Gl0(t)  =  sa9< 


SECTION  C.4.c 
JOINT  TORQUES 

These  torques  are  caused  by  joint  elasticity,  coulomb  friction, 
and  angular  restraints.  The  latter  two  only  arise  if  the  relative 
velocity  of  the  links  about  a  joint  exceeds  some  limit  £,  usually 
about  .01  radians  per  second.  There  are  no  forces  due  to  joint 
torques  in  the  horizontal  and  vertical  direction.  The  equations  for 
the  generalized  forces  in  the  directions  are  as  follows. 

i=9,10  J±(t)  =  0, 

i=8  J8(t)  -  K?[a?(t)  -  otg (t)  +  ag(0)  -  a?(0)]  -J?(t), 

i-4,6  J^t)  =  Ki[ai  ^(t)  -  a±(t)  +  <^(0)  -  <*^(0)]  + 
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±=1,2,5  J±(t)  =  (K^  +  Ki+1)ta1(°)  "  +  Ki+il“i+l(t)-oli+l(0)] 

+  VCtm(t)-ain(0)]  +  J±(t)  “  J±+l(t)’ 

i=7  J7(t)  =  Ky[a8(t)  -  ay (t)  +  a?(0)  -  ag(0)]  -  K1[a7(t)  -a^t)  +  0^(0) 

-a  (0)]  ■!•  Jy(t)  -  Jjtt), 

i=3  J3(t)  =  (K3  +  Ka  +  K5)[o3(0)  -  a3(t)]  +  K3[a2(t)  -  <*2(0)] 

+  KA[aA(t)  -  01^(0)]  +  K^Ia^Ct)  -  0^(0)] 

+  J3(t)  -  J4(t)  -  J5(t), 

where  for  i  =  1-7: 

J1(t)  =  0  for  1^(0  -  ain(t)|<?±, 

J_,(t)  =  C'  sgn[a  (t)-a.(t)]  +  T  (t)  otherwise, 
i  1  m  1  i 

The  specific  form  of  the  angular  restraint  torques  may  vary, 
but  its  general  form  is  shown  in  Fig.  9.  The  simplest  and  most 
usually  used  form  is  linear.  An  actual  restraint  curve  for  the  upward 
motion  of  the  lower  leg  about  the  knee  joint,  as  found  by  J.  W.  Smith 
(ref.  4),  is  also  shown  in  Fig.  9.  A  cubic  fits  this  curve  quite  well. 

Table  6  shows  the  joint  angle  subscript  correspondence. 


i 

1 

5 

7 

2, 3, 4, 6 

m 

- 

7 

3 

8 

i-1 

TABLE  6 


One  of  the  injury  criteria  used  as  an  end— of— program  test  can  be 
applied  here.  This  is  joint  dislocation  due  to  excessive  relative 
angle.  Its  flow  diagram  is  shown  in  Fig.  10. 


35 


37 


SECTION  C.A.d 


1 

’  i 


'X 

SEAT  FORCES 


The  present  representation  of  the  seat  is  that  it  is  composed  of 
a  linear  spring  whose  force  acts  vertically  on  the  upper  or  lower  leg 
at  the  front  edge  of  the  seat,  plus  a  nonlinear  spring  and  linear 
damper  whose  combined  force  acts  vertically  on  the  hip  joint.  This  is 
the  interior  vertical,  not  the  true  direction.  There  is  also  friction 
proportional  to  their  sun  which  acts  horizontally  (also  noninertially) 
at  the  hip  joint.  Logic  in  the  simulation  allows  for  falling  off  the 
seat,  differentiates  between  adults  and  children,  and  prevents  the 
overlapping  of  the  two  types  of  forces.  This  last  is  accomplished  by 
assuming  that,  at  some  distance  z  from  the  front  edge  of  the  seat,  its 
representation  changes  from  this  combination  to  a  linear  spring  whose 
spring  constant  is  dependent  on  the  conditions  at  the  previous  time 
instant . 

Figure  11  depicts  these  various  cases;  figure  12  shows  the  flow 

diagram  by  which  the  net  forces  on  the  body  are  calculated  for  each 

case.  The  parameter  y  has  already  been  defined  and  calculated  in 

zo 

Section  C.2.  The  generalized  force  equations  in  terms  of  these  net 
seat  forces  plus  the  friction  are  as  follows: 

S± ( t )  =  0  for  i  =  1-6, 


8y: 

la. 

ay, 


s7(t)  •  s8(t)  ■  ¥'e ">&• 

/  o 


S  (t)  =  F'(t)-~+  [F'(t)+F  (t)]sin$(t  )  -  f  (t)cos$(t  ), 
7  scix-ss  m  m 
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9y 

S10(t)  =  Fs(t)9yJ+  lF8(t)+F8<t)]cos$<tm>  +  f(t)sln$(tm), 

where,  for  z(t)  j<  L7cos07(t) , 

3yz  2  3yz 

■5-=  z(t)sec  07(t),  0 

7  o 


3yz 

-cos$(t  )tan0_(t),  t —  =  sin$(t  )tan0-,(t) 
m  7  <3y,  m  / 


and  for  z(t)  L7cos  7(t), 

3yz 

=  L7[cos07(t)+sin07(t)tan0g(t)] , 


3yz  2 

=  [z(t)-L7cos07(t)]sec  0g(t), 

8 


3y  3y 

=  -cos$(tm)tan08(t),  ^  =  sin$(tm)tan08(t), 

where  0,  (t)  =  a,.(t)  +  $(t  )• 

11  m 

Another  method  of  preventing  overlap  is  to  change  the  repre¬ 
sentation  of  the  seat  to  be  either  that  of  a  non-rigid  beam  or  that 
of  a  rigid  beam  supported  at  fixed  points  by  the  present  spring  and 
spring-damper  systems.  In  both  cases  the  equations  of  motion  will 
be  different. 

Another  possible  representation  for  the  seat  cushion  is  that  of 
a  non-rigid  beam.  In  particular  it  will  be  assumed  to  consist  of  a 
pair  of  uniform  beams;  one  of  length  l  and  supported  at  both  ends  and 
the  other  of  length  kJ£  cantilevered  at  a  fixed  distance  equal  to  its 
length  from  the  front  edge  of  the  other  beam  and  thus  of  the  seat. 

The  two  beams  and  the  resulting  system  are  depicted  in  figure  13. 
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Figure  13 


The  flow  diagram  is  shown  in  figure  14,  and  auxiliary  variables  and 
parameters  are  calculated  from  the  following  equations: 


z(t)  -  zq  -  x1(t), 

a  -  i-  z(t), 

3  2  2  3  3 

b  =  zJ(t)  -  3k Vz(t)  +  2kV, 

c  =  k£  -  z(t) , 

d  -  ygC  -  yx(t)  +  [^(t)-y^c]  tan  Y0, 


-1  rl"r7 
=  sin 

L7 


.  -l/7  ra. 

=  sin  (T7~)’ 

o 


=  ^  +  d. 


3 

e_  l  6.  e06  A  e_b6 

_  _1 _ A_  *  _  2  z  =  1  z 

^  "  3  3  9  bs  3  3»  *  333 

a  z  (t)  k  Jl  2k  a  z  (t) 


•  Fs  = 


2e0k5,6 . 
2  A 


Body  parameters  involved  are  the  hip  radius  r^,  knee  radius  r^, 
and  backward  ankle  radius  r  .  Beam  parameters,  besides  the  lengths 
mentioned,  are  e^  and  e^,  where  e  is  twice  the  product  of  the  moment 
of  inertia  and  the  Young’s  modulus  for  the  beam. 

Energy  calculations  must  be  done  incrementally  as  for  the  belt 


and  contact  forces. 


There  are  two  reasons  why  this  representation  has  not  been  in¬ 
corporated  into  the  simulation.  The  first  is  that,  since  the  seat 
forces  are  no  longer  derived  from  the  potential  energy,  the  body 


Figure  14 

NON-RIGID  BEAM  SEAT  FORCES  FLOW  DIAGRAM 


equations  would  have  to  be  re-derived.  The  second  is  that  experimental 
confirmation  and  numerical  data  are  as  yet  lacking. 

A  second  possible  alternate  representation  for  the  seat  cushion  is 
that  of  a  rigid  beam  with  some  mass  and  inertia  supported  by  a  linear 
spring  at  the  front  edge  and  by  a  nonlinear  spring  plus  damper  some¬ 
where  near  the  other  end.  Figure  15  shows  the  geometry  and  relevant 
parameters.  It  is  assumed  that  the  horizontal  location  of  the  seat 
center  of  gravity  remains  fixed  at  about  half  its  projected  length. 
Additional  parameters  needed  in  this  case  are  the  unloaded  seat 
cushion  angle  (y  )  and  the  unloaded  y-coordinate  of  the  seat  center 

of  gravity  (y  ),  as  well  as  its  mass  (m  ),  moment  of  inertia  (I), 
gu  s  s 

and  the  distances  of  the  center  of  mass  and  spring-damper  support  from 
the  front  edge  (a  and  a  respectively) . 

The  representation,  however,  requires  two  additional  degrees  of 
freedom  -  those  of  the  y-coordinates  of  the  seat  center  of  gravity 
and  of  its  angle  to  the  horizontal.  These  two  accelerations  can  be 
found  in  terms  of  the  seat  support  forces,  some  parameters,  and  z(t) 
and  y^(t),  so  that  they  are  at  least  uncoupled  from  the  matrix  equation 
and  may  be  solved  for  independently. 

The  equations  are: 

msyg(t)  =  Fs(t)  +  Fs(t)  "  Wo  -  Fso  -  ms8  +  Vl(t) 

I  Y (t)  =  a'[F'(t)-F'  ]  +  (a'-a)F  (t)  +  W  [z(t)-a] 
s  s  so  s  o 

Fs(t)  =  ^emAs(t)  +CsA8(t> 
m=l 
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so 
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RIGID  BEAM  SEAT  GEOMETRY 


Fg(t)  -  sAz(t) 

Ag(t)  =  ygu  -  yg(t)  +  (a-a')[tany(t)-tanYu] 

A  (O  =  y„„  “  y„(t)  +  a’ttany  -tany(t)] 

Z  gu  g  u 

A  (t)  =  (a-a’)Y(t)sec2y(t)  -  y  (t) 
s  g 

The  first  two  determine  the  accelerations;  the  last  five  the  seat 
support  forces  (F  and  F?)  in  terms  of  y  (t),  y  (t),  y(t),  y(t)  and 

S  S  go 

the  parameters  of  the  seat  cg,  s;  a,  aT,  ygu>  Yu) •  At  time  t=0 

there  are  six  initial  conditions  (Fgo,  F*0»  AgQ,  Azo,  yq»  ygQ) 

as  well  as  z  ,  and  parameters  (w  ,  F*  ,  Is,  m  ),  while  all  velocities 

and  accelerations  are  zero.  Thus  there  are  six  operative  equations 

(all  except  the  last,  which  is  satisfied  trivially).  Usually  zq,  Ig, 

m  .  v  •  Y  •  s.  c  .  and  the  3  can  be  measured  and  then  any  four  of 
s  gu  u  s  m 

the  remaining  ten  (a  a',  wo,  F^,  Fgo>  F^,  Ag<)f  A^,  ygQ,  YQ)  will 
determine  the  other  six. 

Figure  16  shows  the  flow  diagram  for  determining  the  forces  on 
the  body  in  terms  of  the  seat  support  forces. 
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Figure  16 

RIGID  BEAM  SEAT  FORCES  FLOW  DIAGRAM 
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SECTION  C.4.e 


BELT  FORCES 


The  restraint  system  consists  of  a  lap  belt  and  a  diagonal  torso 
restraint.  Logic  in  the  program  allows  the  simulation  of  both  belts, 
of  neither  belt,  or  of  either  one  of  them  singly. 

Each  belt's  nominal  deformation  and  rate  is  computed  from  its 
geometry  and  the  present  state  of  the  system.  But  the  nominal 
deformation  is  the  sum  of  the  actual  belt  elongation  and  the  body 
deflection.  For  the  lap  belt  the  latter  occurs  in  the  abdomen;  for 
the  torso  restraint  in  the  chest.  Thus  the  total  deformation  must  be 
partitioned  between  the  belt  elongation  and  the  body  deflection, 
and  likewise  for  the  rate.  Once  the  actual  belt  elongations  and 
rates  have  been  obtained,  they  may  be  used  to  find  the  belt  load 
forces,  which  in  turn  determine  the  generalized  forces.  Recurrent 
loading  with  additive  residual  elongation  is  allowed.  Figure  17 
shows  the  lap  belt  geometry.  The  geometric  equations  for  the 
nominal  deformation  and  rate  are,  for  the  lap  belt  (J-l): 

-1  _ 

♦,(t)  -  sin  iii  j)  i  +  tan 

1  /u^(t)+Vj(t) 

SjXt)  -  /ij(t)  ♦  vj(t)  -  hj‘  -  *l0  +  h1t<t1(t)-e1(t)-4F]» 


8  (t)  -  x^tjcos^t)  +  y^(t)sin^1(t)  -  h^e^Ct), 
where  ^(O-x^O+Aj  and  v^O-y^O+A^  For  the  torso  restraint 
(J»2,3),  the  equations  are: 
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Figure  17 
LAP  BELT  GEOMETRY 
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yo  -  /uju)  +  vj(t)'  -  tJo, 

yt)  »  Xj(t)cosyt)  +  yOiinyt)  +  y  (t)»in[yt)-yt)]; 

where 

^(t)  -  x2(t)  +  h2  C08®2^t^  ~  X3» 
v2(t)  *  y2(t)  +  h2«in02(t)  -  X4, 
u3 ( t)  -  x  <t>  +  h3co.e3(t)  -  X5> 
v3 (fc )  "  y3<t>  +  h3«in63(t)  -  A^. 

The  x^(t),  y^(t)  were  calculated  at  the  beginning  of  section  C.4  and 
the  X ^  in  section  C.2. 

Details  of  recurrent  loading  are  shown  in  figure  18  and  the  flow 
diagram  for  the  belt  forces  in  figure  19,  There  are  essentially  two 
types  of  recurrent  loading,  one  in  which  the  load  is  completely  removed 
before  the  new  loading  cycle  begins  and  the  other  in  which  it  is  not. 

The  former  is  marked  case  1  in  figure  18  and  the  latter  case  2,  In 
both  cases  the  new  cycle  begins  at  the  point  marked  x,  where  5(t) 
equals  zero,  going  from  negative  to  positive.  The  residual  or 
permanent  deformation  for  the  new  cycle  (cn+})  1®  calculated  differently 
in  the  two  cases. 
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case  1. 


Figure  18 

RECURRENT  LOADING  CURVES 
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be  possible  to  neglect  them  and  either  average  the  two  sides  or  use 
one  of  them  consistently. 
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SECTION  C.4.f 


CONTACT  FORCES 

The  contact  surfaces  considered  are  both  sides  of  the  head  rest 
and  seat  back,  the  roof,  windshield,  sun  visor,  upper  and  lower 
instrument  panels,  floor  and  toeboard,  and  steering  wheel  rim  and 
column.  All  of  these  except  the  head  rest,  sun  visor,  and  steering 
wheel  are  assumed  to  be  rigidly  fixed  relative  to  the  interior  of 
the  vehicle.  It  is  further  assumed  that  the  characteristics  of  the 
motion  and  material  of  the  nonrigid  surfaces  can  be  lumped,  with  the 
exception  of  the  column  which  may  be  padded.  In  this  latter  case, 
it  is  assumed  that  the  column  does  not  move  until  its  padding  (if 
any)  has  bottomed  out.  These  assumptions  simplify  the  calculations 
a  great  deal  and  are  sufficiently  realistic. 

In  this  present  simulation,  the  nonrigid  contact  surfaces  are 
characterized  by  load-deflection  relations  rather  than  by  explicit 
differential  equations.  The  latter  is  a  better  and  more  accurate 
representation,  but  requires  that  the  number  of  degrees  of  freedom 
of  the  system  be  raised  by  approximately  six,  since  the  rim  and 
column  would  need  at  least  two  degrees  of  freedom  each  and  the  head 
rest  and  sun  visor  one.  This  will  be  done  in  phase  three. 

Figure  20  shows  the  overall  contact  force  flow  diagram;  its 
more  complicated  parts  are  defined  in  figures  22  to  24  and  in 
tables  7  to  11. 

The  body  reference  point  coordinates  and  their  velocities  are 
calculated  from  the  following  equations ,  referred  to  in  figure  20  as 


equations  f.l* 

i-1,5,  *!(t)  -  x  (t)  x! (t)  -  x  (t) 

6  7:  1  1  1  1 

y[(t)  -  yt(t)  y'l(t)  -  yt(t) 

im2,3:  x|(t)-X3(t)+p|cos[03(t)+Y1]  x£(t)-z3(t)-p.[d3(t)sln[e3(t)+Y1] 

yi(t)-y3(t)+p^8in[03(t)+T1]  y’ (t)-y3(t)+p|03(t)cos[03(t)-hr1] 

1-4:  x£(t)  -  x^(t)+p^cos04(t)  x^(t)-x^(t)-p^0^(t)«ln04(t) 

y4(t)  “  y4(t)+P4sln04(t)  y^(t)-y4(t)+p404(t)cos@4(t) 

1*8:  xg(t)“x^(t)+LgCOS0g(t)  xg(t)-x7(t)-L8Sg(t)sln08(t) 

y8(t)-y7(t)+L8slneg(t)  y^(t)-y7(t)+Lg0g(t)cos08(t) 

1-9:  x’(t)  -  x?(t)  i’(t)  -  x?(t) 

y’<t)  -  y?(t)  y*(t)  -  y?(t) 

l-iO:  x'  (t)-x  (t)+L  cos0  (t)  x'  (t)-x  (t)-L  0  (t)sln0  (t) 

10  666  10  666  6 

yio(t)“y6(t)+L6Sine6(t)  yjQ(t)-yg(t)+L^0^(t)cos0^(t) 
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Table  7  shows  what  sur fsccr;  cav  hr-  contacted  by  which  body 
segments  for  the  three  occupant.  poriti.o.r 


Body 

Segment 

Subscript 

i 

Driver 

hip 

1 

> 

back 

2 

9 

chest 

3 

8,10 

head 

4 

1 , 2 . 3,8,9 ,10 ,11 

shoulder 

3 

10 

elbow 

6 

9 

knee 

7 

5 

foot 

8 

6,12 

shin 

9 

hand 

10 

l_J‘c  l  _ _ 

sub script  k 

Fiiov’l  -Scat 
?fr,r>rrp;er 

9 

o 

1.2,  V  ,9,11 
9 

K 

5,6,12 

f>> 

/t.b  J1 


Rear-Seat 

Passenger 

9* 

9* 

7 

1,7,9*, 11,1* 

9* 

7 

6*, 7,12* 

_liii _ 


where  1*  indicates  the  l  of  the  fro.'!;  t:  f-  head  rest,  5*  the  lower 


edge  of  the  lower  panel ,6*  the  r 'loping  r cation  of  the  rear-seat  floor, 
9*  the  rear- seat  back,  10  *  the  eld*.  rin*  of  the  steering  wheel,  and 


U*U: 


/  c 


12*  the  rear  seat:  floor 


Equations  f.2  for  T  <t),  6lk(t)  and  6lk<t),  as  derived  from 
the  various  geometries  shown  in  figure  21,  are  as  follows: 

6ik(t)  -  rt  +  [x[<t)-x£]sini|ik  +  [yk-y[(t))cos*k, 

«lk(t)  •  y^(t)cos<|<k  -  x|(t)sin4ik, 

Tt(t)  -  [x^-xj(t)]cosi(ik  +  [yk-yj(t)]sin^k, 
for  the  (i,k)  sets  of  table  8,  with  T^t)  omitted  for  k-11,12.  Special 
relations  needed  here  are:  xJVxjJ-D^cosil^*  y4“y2“D28in^2’  ,  ♦jj"180# , 

y'll“y2’  ''’ll"0''  yi2"y6“D68in,|,65 

6  (t)  -  +  r.  -  [x* (t)-x"]cosY  +  [y' (t)-y"]sinY, 

ik  21441  41 

«lk(t) '■  x^(t)c°SY  -  y^(t)sinY, 
for  i«4,  k«l  and  where 

y  ■  tan 

6lk(t)  ■  r7  -  [x^-DjCosi(»j-x^ (t)  ] sin [ ®g (t)-oQ] 

+  [yj-D5sini|<5-y’(t)]cos[e8(t)-oo] , 

«ik(t)  -  y2(t)cos[68(t)-ao]  -  x^(t)sin[08(t)-ao]  +  T9(t)eg(t), 

Tt(t)  -  [x^-D5cosi)<5-x^(t)]cos[08(t)-aol-[y^-D5siniJ>5-y^(t)]sin[08(t)-aoJ> 

for  i*9.  k«5*.  where  a  was  defined  and  calculated  in  section  C.2. 

o 
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subscript 

subscript  k 

i 

Driver 

Front  Pass. 

Rear  Pass. 

1 

9 

9 

9* 

2 

9 

9 

9* 

3 

8,10 

4 

7 

4 

2,3,8,9,10, 

11  2,3,4,9,11 

1*»7,9*,11 

5 

10* 

- 

- 

6 

9 

9 

9* 

7 

5 

5 

7 

8 

6,12 

5,6,12 

7,12* 

10 

5.10*.ll 

_ *,5,11 _ 

Ml _ 

TABLE  8. 


Figure  22  defines  the  T^t)  contact  check  and  also  takes  care 
of  those  cases  for  which  T^(t)  is  omitted. 

Table  9  specifies  those  combinations  of  i  and  k  which  are  termed 
undesirable f  for  use  in  the  undesired  contact  end  of  program  test. 


subscript 

i 

Driver 

subscript  k 

Front  Pass.  Rear  Pass. 

3 

8,10 

4 

7 

4 

2,3,8,10,11 

2,3,4,11 

1*,7,11 

5 

10* 

- 

- 

7 

5 

5 

7 

8 

- 

5 

7 

9 

5* 

5* 

- 

10 

,  .  5, IP.*., .1.1 _ 

4,5,11.  _ 

7.11 _ 

TABLE  9. 
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is  k 
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(t)  CONTACT  CHECK  FLOW  DIAGRAM 


Table  10  defines,  both  numerically  and  descriptively,  the 
composition  of  any  particular  contact  surface,  whether  a  single 
material,  a  padded  surface,  or  a  nonrigid  surface. 


Contact  Surface 

k 

i 

Composition 

head  rest 

1 

8 

non-rigid  motion 

windshield 

2 

2 

glass 

sun  visor 

3 

9 

non-rigid  motion 

upper  panel 

4 

3,6 

padded  metal 

lower  panel 

5 

3,6 

padded  metal 

toeboard 

6 

4,5 

mat  on  floor 

back  of  front  seat 

7 

1 

seat  material 

steering  column 

8 

3,10 

padding  plus  non-rigid  mot 

seat  back 

9 

1 

seat  material 

steering  wheel  rim 

10 

7 

non-rigid  motion 

roof 

11 

3,6 

padded  metal 

floor 

12 

4^5_ 

mat  on  floor 

Body  Segment 

i 

b 

for  non-rigid  motion 

abdomen 

1 

11 

chest 

3 

12 

face 

4 

13 

shoulder 

5 

14 

TABLE  10. 

Table  11  shows  the  combinations  of  i  and  k  for  which  contact 
deflection  is  shared;  i.e.  deflection  occurs  both  in  the  body  and  in 
the  surface. 


subscript 

i 

Driver 

subscript  k 
Front  Pass. 

Rear  Pass. 

3 

8,10 

4 

7 

4 

8,10 

4 

1*,7 

5 

10* 

- 

- 

TABLE  11 
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Out 


Figure  23 

SHARED  DEFLECTION  CONTACT  FORCE  FLOW  DIAGRAM 
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Figure  23  displays  the  flow  diagram  for  this  shared  deflection 
force,  while  figure  24  contains  that  for  deflection  in  or  of  the 
surface  only.  Figure  25  depicts  a  bi-material  force-deflection  curve. 

Equations  f.3  define  the  tangential  velocity  of  the  body  along 
the  contact  surface,  as  follows: 

v  m  x'  (t)cost|i  +  y'(t)sin4  +  r  0  (t)  for  the  Table  8  (i,k), 

T  i  k  i  ^  xi 

*  y^(t)cosy  -  x^(t)siny  +  ^or  k*l, 

VT  “  iy(t)cos[08(t)-oo]  +  y^(t)sin[68(t)-ao]  + 
for  i»9, 

Equations  f.4  define  the  present  horizontal  and  vertical  contact 
force  components ,  relative  to  the  interior  of  the  vehicle,  for  each 
contacting  body  segment  in  terms  of  the  normal  and  tangential  forces 
and  the  angle  of  the  contact  surface  with  respect  to  the  floor.  The 
program  also  allows  for  the  case  in  which  one  body  segment  contacts 
more  than  one  surface  at  the  same  time,  which  is  theoretically 
possible  but  practically  unlikely.  Special  relations  needed  in  these 
equations  are:  ^-180°,  <|/12-0°.  The  proper  (i,k)  sets  from 

Table  7  must  be  used. 

?S(t)  -  P  (Osimli,  +  P'(t)cosi|> 
i  i  k  l  K 

PC^t)  -  P1(t)co8i|!k  -  P^(t)sini|ik, 
i-9:  PS9(t)  -  P9(t)sin[ao-68(t)]  +  P£(t)cos [ao-08(t) ] 

PC9(t)  -  P^(t)sinl08(t)-ao]  -  P9(t)cos[eg(t)-ao]. 
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After  all  these  contact  force  components  have  been  calculated 
for  all  the  contacting  body  segments,  the  three  main  force-based 
injury  criteria  can  be  evaluated.  Figure  26  shows  the  flow  diagram. 

These  are  the  tension  dislocations  of  joints  and  compression  and  shear 
fractures  of  bones.  Table  5  in  section  C,1  shows  the  correspondence 
between  the  subscript  i  and  its  particular  bone  in  each  case.  The 
forces  are  calculated  in  equations  f.5,  f . 6 ,  and  f.7  respectively. 

Equations  f,5 

FT  -  (F  +PC  +PC  +PC  -F'-PC  -PC  )sin0 

1  s  1  2  4  s  8  9  i 

+  (f-PS1-PS2-PSA-Ft+P37+PSg+PS9)cosei, 

FT  -  (PC  +PC  -F  )sin6  +  (PS  -PS  -PS  -f )cos0 

2  24s  2  1^4  2 

+  F^cos  ($2^3)  ~  F2cos(e2-tJ>2)  , 

FT  -  (PC  -F  -PC  ) sin6  4*  (PS  -f-FS  )cos0 

3  4  s  1  3  1  /  3 

+  F2cos(Q3~*2)  +  F3cos  (63-4)3)  , 

FT  -  (PC  -F  -PC  -PC  ) sin 0  +  (PS.+PS  -f-PS . )cos0. 

4  4sl2  4  12  4  4 

+  F2cos(e^~4>2)  -  F1cos(e4-^1)  -  F3cos(eA-4)3) , 

FT  -  (PCc+PC  +PC  -F  -PC.-PC  -PCjsine^  -  F.  cos  (6,-4*!  ) 

5  5  6  10  s  1  2  4  5  i  ->  1 

+  (PS1+PS2+PS4-f-PS5-PS6-PS10)cos65 

+  F2cos(65*-<^2)  -  F3coc(03-4>3)  , 

FT6  «  (PC6+PC10-PC2-PC5)sin66  +  (PS2+PS5-PS6-PS1())  cos06 , 

FT?  -  (PC7+PC8+PCg-Fs-F^)sine8  +  (PS^f-F^-PSj-PSg-PS^cosOg. 
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FORCE- INJURY  CRITERIA  FLOW  DIAGRAM 


Equations  f.6 


FCj  “  (f-PS1+PS2+PS4)cos01  +  (F8+PC1-PC2-PC4)sinel 
+  F2C08  (e^j)  -  F3C08(61-4i3)  , 

FC2  -  (f-PS1+PS2+PS4)cos02  +  (F8+PC1-PC2-PC4)sin02 
+  F1cos(02-01)  +  F3cos(02-<|>3), 

FC3  -  (f-PS1+PS4)cos03  +  (F8+PC1-PC4)sin03  +  F^os  <0  -@3) 

-  F2cos(03-(|>2), 

FC^  omitted, 

FC5  -  (PS6+PSlo-PS2-PS5)cos05  +  (PC2+PC5-PC6-PClo)sin05, 

FC6  "  (PSi0"PS5-PS6)cOse6  +  (PC5+PC6-PC10)8lne6* 

FC?  -  (f-PS1+PS7+PSg+PSg)cos07+  (Fs-PC7-PC8-PCg)sin67, 


FC 

8 


(f-PS1-PS7+PS8)cos0g  +  (Fs+F^+PC7-PC8)sin08. 


Equations  f.7 


FS. 


FS. 


FS. 


FS 


FS. 


FS. 


FS. 


FS. 


FS. 


F2sin(02-<|>2), 
F3sin(03-()>3) , 

V 

-PS2sin03  -  PC2cos03, 

-PS^sinQ^  -  PC3cos63, 

F'cosQ.,, 
s  7 
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Finally  the  contact  generalized  forces  are  calculated  from 
Equations  f.8,  as  follows: 

6,10  6,10 

Q. (t)  -  L  [sin0, (t)  l  PS  (t)+cos0  (t)  £  PC  (t)]  -  r  Pj(t), 

1  1  1  i-1  1  i-1 

6,10  6,10 

Q-(t)  ■  L  [sin0  (t)  l  PS  (t)+cos0  (t)  £  PC  (t)], 

2  2  L  i«*2  1  i»2  1 


Q^(t)  -  L3[PS4(t)sin93(t)+PC4(t)cose3(t)]  -  r2P^(t)  -  r3P^(t) 


6,10  6,10 

+  L.[sin0o(t)  >  PS  (t)+cos8  (t)  2,  (t)] 

4  3  i-5  1  3  i-5  1 

3 

+  l  p^{PS1(t)sin[0i(t)+Yi]+PC1(t)cos[63(t)-hri]}, 

Q  (t)  -  p^[PS4(t)sin04(t)+PC4(t)cos04(t)]  -  r4P^(t), 

6,10  6,10 

Q  (t)  -  L  [sine  (t)  l  PSi(t)+cos05(t)  l  PC^t)]  -  ^’(t)  -  rfiP£(t>. 
3  D  i»5  i®5 

Q6(t)  -  L6[PSlo(t)sin06(t)+PClo(t)cos06(t)]  -  riopio(t)’ 

9  ? 

Q7(t)  -  L?[sine7(t)  l  PSi(t)+cos07(t)  l  PCi(t)]  -  r?P^(t), 
i«7  i*3? 

Qg(t)  “  L8[PS8(t)sine8(t)+PC8(t)cos08(t)]-r8P^(t)-r7P^(t)-T9(t)P9(t). 
10 

Q  (t)  “  l  [PC  (t)sin*(t  )  -  PS  (t)cos*(t  )], 
y  i»l  1  m  x 


Q  (t)  -  l  [PS  (t)»in*(t  )  +  PC  (t)cos$(t  )]* 
10  i«l  1  m 
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SECTION  C . 5 


ACCELERATIONS 

In  the  present  simulation,  the  vehicle  accelerations  are  read 
from  tables.  These  include  the  horizontal  and  vertical  accelerations 
of  the  vehicle  center  of  gravity  and  its  pitch  acceleration. 

The  accelerations  of  the  generalized  coordinates  for  the  body 

-V  - V 

are  found  by  inverting  the  matrix  equation  Az  =  b  and  solving.  The 
vector  £  has  already  been  computed  in  Section  C.4.  The  matrix  A 
is  shown  in  figure  27,  where  the  a  were  defined  and  computed  in 
Section  C.2. 
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and  of  the  individual  body  segments: 

KE±(t)  =-|m.[uJ(t)-K^(t)]  +^1^(0 

for  i  -  1  to  8,  where 

Ut(t)  =  ih(t)  -  I  ^(OrillB^t)  - 

m-1 


i=l-4 


i-l 


Ut)  =  yh(t)  +  I  Vm(t)cOSam(t)  H 

m=l 


2 

/U  (t)  =  ^(t)  -  l  Lmam(t)sinam(t) 
{  m=l 


i-l  . 

-  T  L  a  (t)sina  (t)  -p.a. 
^  m  m  m  l  i 

m=5 

V±(t)  =  yh(t)  +  Iu(t)c°3V(t) 


i-l  . 

+  T  L  a  (t)cosa  (t)  +  p.a 
Lc  m  m  m  i  i 

m=5 


i-l 

u.(t)  =  ih(t)  -  l  Lm;m(t)sinam(t) 

x  m=7 


i=7 ,8 


i-l 

V  (t)  =  yh(t)  +  I  Lman.(t)cosam(t) 
m=7 


Pia±(t)sinai(t) , 


p .a. (t)cosa .  (t) , 
li  i 


-  L^a^(t)sina^(t) 
(t)sinai(t) , 


+  L4a3(t)cosa^(t) 
(t)cosa^(t) , 


-  p.a.  (t)sina  (t) 
li  i 


+  p  .a. (t)cosa . (t) 
ii  i 
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The  conserved  and  dissipated  energies  in  the  seat  are  due  to 
its  springs,  linear  and  nonlinear,  and  to  the  damper  and  friction; 
i.e. : 


„  mi*  1  /  .\ 
(t) 


CE(t)  =  \  s[y1(t)  +  yz(t)-2otanYo]2  +  \  I  " “hi 

m=l 


DE(t)  =  DE(t-At)  + -|{cs[y1(t)+y1(t-6t)][y1(t)-y1(t-At)] 

+  [f (t)+f (t— At) ]  [x^ (t)-x^(t-At) ] } . 

In  the  joints,  the  conserved  energy  is  due  to  the  elasticity 
and  the  dissipated  to  friction  and  restraints;  thus  for  each  joint 
(i  =  1  to  7): 

CEi(t)  =  K.[0m(t)-0i(t)-H3.o-0mo]2, 


DE.(t) 


IVtMi(t>  +  8lo  -  °J  +  2 J0tTi(C> 

+T . (t-At) ] [0  (t)-0  (t-At)-0  .  (t)+0, (t-At) ] 

l  mm  ii 

for  |0m(t)  - 


0  o 


therwise. 


The  gravitational  energy  is  all  conserved;  for  each  body  segment 
(i=  1  to  8): 

CEi(t)  =  gmi[y^g(t)cos$(tm)  -  y^cos<i> (t) ] . 

where 

y^8 (t)  «  y±(t)  +  Pisin0i(t)  for  i  =  1  to  6, 
yyg(t)  *  y±(t)  +  P7sin07(t), 
ygS(t)  =  y7(t)  +  pgSin0g(t) . 
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The  conserved  energy  for  each  contact  (i  =  1  to  10)  is  CE^Ct) 
R  E.  (t)  (for  properil);  the  dissipated  energy  is 

AS  1 

DE1(t)  =  (l-R^)E±(t) . 

The  belt  conserved  energy  Is  (for  j  =  1  to  3) : 

CE^t)  =  Ecj(t), 
dissipated  Is: 

/i-R .  •> 

DEj  (t)  = -L^“LL-Ecj(t). 

The  last  piece  of  output  data  at  present  is  the  contact  force 
position,  defined  by  T^t)  for  i  =  1  to  8. 
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SECTION  C.  7 


END-OF-PROGRAM-TEST 

The  end-of-program  option  allows  the  selection  of  several  tests. 
At  present  these  choices  are  time  only,  time  plus  undesired  contact, 
and  time  plus  injury.  Thus  the  time  duration  of  the  program  run  is 
always  one  of  the  criteria  for  the  end-of-program  test,  and  in  fact 
is  the  only  one  contained  in  this  section.  It  is  shown  in  figure  2ft. 
The  undesired  contact  test  is  contained  in  Section  C.4.f  as  Table 

9  and  in  figure  20. 

The  injury  test  occurs  in  several  places;  that  for  joint  dis¬ 
location  due  to  excessive  relative  angle  in  Section  C.4.c  in  figure 

10  and  the  others  in  Section  C.4.f,  figure  20  and  figure  26. 

Another  criterion  that  could  be  used  as  an  end-of-program  test 
is  the  cessation  of  body  motion. 

Which  end-of-program  test  is  chosen  would  depend  in  large  part 
on  what  the  simulation  is  being  used  for.  If  the  purpose  is  to  in¬ 
vestigate  or  test  various  restraint  systems,  then  time  plus  undesired 
contact  would  probably  be  the  best  criterion  to  use.  If  the  interior 
vehicle  design  is  the  object,  then  the  time  plus  injury  criteria  would 
be  of  value.  If  the  purpose  is  just  to  find  the  motion,  then  time 
alone  will  be  sufficient. 
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Figure  28 

TIME  END-OF-PROGRAM  TEST 


SECTION  8 


MODIFICATIONS 


It  would  be  extremely  useful  for  the  program  itself  to  be  able 
to  evaluate  its  own  results,  change  certain  input  parameters  or  data, 
and  then  initiate  a  new  run.  The  criteria  for  such  evaluations  and 
changes  would  depend  on  the  ultimate  use  of  the  program,  as  would  the 
selection  of  parameters  to  be  varied. 

For  example,  if  the  program  were  being  used  to  evaluate  restraint 
systems,  then  the  natural  parameters  to  be  varied  would  be  those  of 
the  belt  material  and  their  anchorage  locations.  The  evaluation  criteria 
could  be  the  use  of  the  undesired  contact  end-of-program  test,  the 
injury  criteria  for  the  torso  and  neck,  and/or  large  forces  in  the  belts 
themselves. 

The  method  of  changing  the  parameters  could  be  either  systematic 
or  random,  and  they  could  be  changed  either  singly  or  all  at  the 
same  time. 
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SECTION  D 


A  SIMPLE  ANALOG  TO  COLLISION  BETWEEN  AUTO  OCCUPANTS 

The  following  describes  a  simple  mathematical  analog  to  the 
collision  between  front  and  rear  auto  occupants.  The  front  seat 
occupant  is  restrained  by  a  seat  belt  and  shoulder  harness  whereas 
the  rear  seat  occupant  is  restrained  only  by  a  seat  belt.  The 
simple  model  is  shown  in  the  following  sketch. 


Points  A  and  B  represent  the  hip  joints  about  which  the  upper  portion  of 
the  bodies,  Wj^  and  w2  pivot  due  to  the  seat  belt  restraint.  The  rear 
seat  occupant  is  restrained  from  pivoting  only  by  the  front  seat  back. 
For  simplification,  this  is  represented  by  a  rigid  stop  when  02  =  60°. 
The  front  seat  occupant  is  restrained  by  a  shoulder  harness  with 
spring  and  damping  constants,  k^  and  c.  The  angle  can  reach  a 
may-iTTnim  of  -  14°.  The  forward  occupant  is  also  restrained  by  a 
hip  joint  torque,  T.  It  has  been  found  that  the  inclusion  of  this 
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restraint  on  the  rear  occupant  increases  the  amount  of  time  between 
initial  impact  and  final  stop  by  about  8%.  The  initial  forward 
velocity  is  v.  The  objective  of  the  following  is  to  show  a  simple 
analog  to  the  difference  between  interoccupant  motion  in  the  cases 
where  the  front  shoulder  harness  is  an  undamped  spring  and  where  it 
is  critically  damped. 

The  equation  of  motion  for  the  front  seat  occupant  is 
I0  G1  +  *1C®1  +  ^l^l^l  ±  T  =  0 

The  plus  sign  on  T  applies  when  0  is  positive.  The  undamped  solution 
for  initial  conditions 

01(O)  =  0,  0X (0)  =  tu 

is  given  by 


For  the  critically  damped  case, 

Sl  c 

e  -  ±  -rr-exp(-  -f-t)  +  <u  *  2k“l-)t  exp(" 

1  i  ^^0  10 

The  equations  of  motion  and  solution  for  the  rear  seat  occupant 
under  the  same  initial  conditions  are 

V2  '  0 

0  =  wt 

To  make  the  situation  as  realistic  as  possible  the  system  parameters 

% 

have  been  determined  as  follows.  The  two  occupants  are  assumed  to  be 
large  (95th  percentile  males)  with  a  weight  of  the  upper  part  of  the 
body  of  about  110  lb  and  a  length  of  about  38  inces.  The  distance  from 


V 

K' 


+  -X- 
Vi 
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the  pivot  point  to  the  center  of  gravity  was  taken  to  be  17  inches 
yielding  1^  =  117  lb- in-sec^. 

The  initial  angular  velocity  was  found  by  noting  that  initially  each 
body  possessed  a  given  linear  momentum  which  had  a  moment  of  momentum 
about  the  pivot  point  of 

mvp  =  IqO ) 

For  a  20-mph  collision,  co  =  14.7  rad/ sec. 

By  trial  and  error  it  was  found  that  a  spring  constant  K  =  25,000 
lb/in.  stopped  motion  of  the  forward  occupant  in  0.026  sec  for  the 

undamped  case  at  about  0^  =  14  .  Using  the  previously  determined 

% 

parameters  the  damping  coefficient  was  found  to  be  c  =  830  lb-sec. 

The  motion  in  the  undamped  case  preceded  as  follows.  At  t  -  0.026 
sec,  the  forward  thrust  of  the  front  occupant  had  been  arrested  by  the 

'Xj  ^ 

shoulder  harness  at  0^  =  14  .  The  rear  occupant  had  pivoted  to  0^  =  22 
in  this  time.  The  complete  time  history  is: 


time (sec) 

(approx) 

©2 (approx) 

0 

0 

0 

0.026 

14° 

22° 

0.050 

1° 

42° 

0.068 

-  8° 

60° 

0.075 

-13° 

60° 

The  simple  sketch  below  shows  that  contact  must  have  occurred  between 
the  two  occupants  between  t  =  0.050  and  t  =  0.075  if  a  typical  seat 
separation  (L  =  34  inches)  is  used. 
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Dotted  lines  represent  the  final  positions. 

For  the  damped  case  numerical  calculations  show  that  the  solution  is 

dominated  by  the  term 

HCt 

0  =  wt  exp  — 

This  term  is  never  negative,  and  thus  the  springback  causing  interoccupant 

Of  0 

collision  using  this  model  is  impossible.  By  the  time  t  =  0.125,  0^-1 
The  major  objective  of  the  simple  study  has  been  accomplished.  The 
model  has  shown  that  collision  between  the  two  masses  may  or  may  not 
occur  depending  on  system  damping.  Other  parameters  in  the  model 
besides  damping  will  also  play  a  major  role.  These  are  size  of  masses, 
pivot  point  location,  rate  of  speed,  geometry  of  restraints,  pivot 
point  restraint  friction,  etc.,  all  of  which  are  in  the  model  and  can 
be  varied.  These  parameters  can  be  considered  analogous  to  size  of 
occupant  (not  all  drivers  and  occupants  are  95th  percentile  male 
dummies),  seat  location,  rate  of  speed,  harness  and  seat  belt  geometry, 
joint  friction,  and  seat  design  and  properties. 
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